Acoustic black holes in a two-dimensional "photon fluid 
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Optical field fluctuations in self-defocusing media can be described in terms of sound waves in a 
2D photon-fluid. It is shown that, while the background fluid couples with the usual flat metric, 
sound-like waves experience an effective curved spacetime determined by the physical properties of 
the flow. In an optical cavity configuration, the background spacetime can be suitably controlled by 
the driving beam allowing the formation of acoustic ergoregions and event horizons. An experiment 
simulating the main features of the rotating black hole geometry is proposed. 
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INTRODUCTION 

Progress in understanding critical phenomena in gen- 
eral relativity and quantum field theory in curved space- 
time, suffers greatly from the lack of experimental feed- 
back. For this reason, strong efforts have been made 
in order to find non-relativistic systems, experimentally 
testable in the laboratory, able to simulate gravitational 
spacetime geometries and cosmological solutions [l| . The 
first proposal dates back to 1981, when Unruh showed 
that the propagation of sound waves in an inhomoge- 
neous flowing fluid is analogous to that of a massless 
scalar field in curved Lorentzian spacetime [2[ . A sur- 
face in this flow, where the normal component of the 
fluid velocity is inward pointing and equal to the local 
speed of sound behaves as the event horizon of a grav- 
itational black hole, i.e. sound waves cannot propagate 
through this surface in the outward direction. This kine- 
matical analogy implies that acoustic black holes should 
emit Hawking radiation 0] in the form of a thermal bath 
of phonons [1, 0] , thus making possible the experimental 
verification of such phenomenon. 

Following Unruh's prediction many different kinds of 
analogue models has been proposed but Bose- 

Einstein condensates (BEC) [1, 0] are currently the best 
workbenches for this kind of investigation. The exper- 
imental implementation of acoustic black holes in such 
systems could be possible by setting suitable external 
potentials to obtain the required background flow. In- 
teresting proposals have been done also in the optical 
domain when Leonhardt and Piwnicki proposed to cre- 
ate an optical black hole The idea is feasible since 
the propagation of light in a moving medium resembles 
many features of a curved spacetime and very recently 
ultrashort pulses in optical fibers demonstrated the for- 
mation of an artificial event horizon 11]. 

It is the aim of this work to show that exploiting 



where no is the linear refractive index, n<i > is the mate- 
rial nonlinear coefficient and \E\ 2 is the optical intensity. 
As a consequence light propagating in a self-defocusing 
medium induces a local negative bending of the refrac- 
tive index which, in turn, affects the light beam itself. 
At a microscopic level this can be described in terms 
of an atom-mediated repulsive interaction between pho- 
tons which leads to the formation of a "photon-fluid" 
[liL Il3j |. It will be shown that linear excitations of 
such fluid (sound waves) propagate in an effective curved 
spacetime determined by the physical properties of the 
optical flow. However, the simple propagation in a self- 
defocusing medium does not allow the implementation 
of the optical background required to generate a station- 
ary "black-hole" geometry. In contrast, if the medium is 
placed in an optical cavity, suitable choice of the injected 
beam profile leads to the formation of acoustic ergore- 
gions and event horizons, making possible the experimen- 
tal simulation of an acoustic rotating black hole. The 
possibility to control the spacetime geometry together 
with the fact that light signals are much easily detectable 
than acoustic disturbances in fluids make the system here 
proposed an interesting alternative for "analog-gravity" 
experiments. 



ACOUSTIC METRIC IN SELF-DEFOCUSING 
MEDIA 

First consider the simple propagation of a monochro- 
matic optical beam of wavelength A, in a self-defocusing 
medium (without cavity). The slowly varying enve- 
lope of the optical field follows the well-known Nonlinear 
Schroedinger Equation (NSE) [3] 
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the relation between nonlinear optics and fluid dynamics 
rotating acoustic black holes can be created in a self- 
defocusing optical cavity. Self-defocusing media possess 
an intensity-dependent refractive index, n = uq — ?i2|-E| 2 , 



where z is the propagation direction, and k = 27rno/A 
is the wave number. W 2 E, defined with respect to the 
transverse coordinates (x, y), accounts for diffraction, 
while the nonlinear term describes the self-defocusing ef- 
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feet. The link with fluid dynamics becomes evident writ- 
ing the complex scalar field in terms of its amplitude and 
phase, E = p 1 l 2 e % $ . This representation converts Eq. (Q]) 
into the hydrodynamic continuity and Euler equation 
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where the optical intensity p corresponds to the fluid den- 
sity, v = j^-V0 = Vt/i is the fluid velocity (here c is 
the speed of light) and the propagation direction acts as 



time t 



Apart the optical coefficients, Eq. (2) 



(3) are identical to those describing the density and the 
phase dynamics of a 2D BEC in the presence of repul- 
sive atomic interactions [lH (the interaction is attrac- 
tive if self- focusing nonlinearity is considered) . The opti- 
cal nonlinearity, corresponding to the atomic interaction, 
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provides the bulk pressure P = c " 2 : f while the last term 
in Eq. (3) (quantum pressure) has no analogy in classi- 
cal fluid mechanics. In optics, it is a direct consequence 
of the wave nature of light (it arises from the diffraction 
term) and is significant in rapidly varying and/or low in- 
tensity regions such as dark solitons core and close to 
boundaries. 

The evolution equation for sound waves in the photon 
fluid can be obtained by linearizing Eqs. (2)-(3) around a 
background state, since acoustic disturbances are defined 
as the first order fluctuations of the quantities describing 
the mean fluid flow. By setting p = po + epi + 0(e 2 ), and 
if) — -00 + e0i + 0(e 2 ) we obtain 
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When the quantum pressure is negligible, Eqs. (4)-(5) 
can be reduced to a single second order equation for the 
phase perturbations 



ih ( ^ (dtVi+Vo-Wi) 



+V • po V^i - 



Pov 



(c^+vo- Wi) =0. (6) 



where c s is the local speed of sound, usually defined as 
Note that Eq. ([6]) has the form 
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of a wave equation for a massless scalar field A-0i = 
7=J^ (V~9 9^ dui>i) = propagating in a (2 + 1)- 
dimensional curved spacetime whose (covariant) metric 
reads 
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where g — detlg^) and I is the two-dimensional identity 
matrix. Equivalently, the line element in polar coordi- 
nates on the plane is 

ds 2 = (—) 2 [-(c 2 s -v 2 )dt 2 -2v r drdt-2v e rd8dt+dr 2 + (rd0) 2 

(8) 

where u 2 . = v 2 + v 2 , v r — dripo and vg = ^detpo- The 
acoustic metric (7), already found in BECs [8[ and clas- 
sical fluids 0, 0|i can have ergoregions and event hori- 
zons depending on the physical properties of the flow. 
The region in which the fluid velocity exceeds the speed 
of sound, i.e. Vq > c 2 defines the ergosphere, where no 
physical objects can remain at rest relative to an inertial 
observer at infinity. The event horizon, instead, is de- 
fined by the surface where the radial component of the 
fluid velocity v r equals c s . Note however that this geo- 
metrical interpretation fails when the quantum pressure 
cannot be neglected, for instance when the length-scale 
of spatial variations of the density is smaller than the 



critical length £ 



This point can be clarified 
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considering regions of nearly homogeneous background, 
where pi and ip\ can be treated as slowly varying am- 
plitude plane waves (eikonal approximation) Then 
Eqs.(4)-(5) lead to a curved-space generalization of the 
Bogoliubov dispersion relation in a Bose gas [l6| 
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where K is the wave number of the sound mode and f2 
its angular frequency in the laboratory frame. Assum- 
ing vq = (no background velocity) it is immediate to 
see that modes of wavelength A ^> £ follow the standard 
phononic dispersion relation, Q w c s K, while at high- 
energies (A <C £) the quantum pressure contribution is 
dominant and fl w 2ET^- 2 - ^ n ^ ms case the group veloc- 
ity can grow without bound allowing short-wavelengths 
modes to escape from behind the horizon. The length £ 
(the analog of the "healing length" of BEC) is the scale of 
the violation of Lorentz invariance, usually identified in 
quantum gravity phenomenology with the Planck scale. 

In the long- wavelength limit, where the picture of an 
effective metric makes sense [f9[ , acoustic black holes can 
form in presence of a suitable optical background which, 
however, is constrained to solve Eqs. (2)-(3). One of 
the most studied fluid flows to pursue "analog gravity" 
scenarios is the stationary vortex flow, being the hy- 
drodynamic structure closer to rotating black holes [20(. 
Its optical counterpart arises from the self-trapping of 
a phase singularity embedded in a broad optical beam 
due to the counterbalanced effects of self-defocusing and 
diffraction [3]. The resulting structure (vortex soliton) 
is characterized by a dark core and a specific helical wave 
front, Eq = p]/ 2 (r)e l ^'°, where %pa — m6 and po(0) = 0. 
However, while such flow geometry possesses an ergosur- 
face when c s (r) = vg = j ^ n r an event horizon cannot 
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form since v r is identically zero. As remarked by Visser 
2l|| , provided that c s remains positive, an hydrodynamic 
event horizon requires the vortex to have a central sink, 
i.e. to possess a non-zero inward radial velocity (collaps- 
ing vortex) [2^ | . Unfortunately collapsing vortices appear 
only as transient solutions of Eq. flj and once the ini- 
tial conditions has been chosen their dynamics depends 
only on the nonlinear wavenumber shift in the medium. 
For this reason, a more controllable system in which the 
photon-fluid can be created and allowing the choice of 
the background flow profile is required. In the next sec- 
tion, it will be shown that both these characteristics are 
fulfilled in a self-defocusing optical cavity. 



ACOUSTIC HORIZONS IN AN OPTICAL 
CAVITY 

In the mean field approximation, the slowly varying 
envelope of the intracavity field follows a damped-driven 
version of the NSE 
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where L is the cavity length, to is the pump frequency, 
6 = lu — Lu cav is the detuning between the pump and the 
cavity resonance, Ed is a coherent driving field, which is 
proportional to the incident field, and T — cT /2tlqL is the 
cavity decay rate, where T is the mirrors transmissivity. 
It is immediate to see that Eq. (TTOl) is a dissipative-forced 
system, where the losses —iTE are balanced by the co- 
herent injected field iTEd- Therefore, if T is finite, the 
presence of the driving field " pins" the intracavity phase 
profile allowing the external control of the photon-fluid 
background, and thus the generation of "black- holes" ge- 
ometries. At the same time, if T is sufficiently low the 
photons remain trapped inside the cavity long enough 
so that a thermalized condition is achieved after many 
photon-photon interactions, thus allowing the formation 
of the photon fluid Q . As a further advantage note that, 
contrary to Eq. (1) where the variable z acts as time, in 
the cavity configuration the field envelope evolves in " real 
time" t, thus allowing experimental studies of real sound 
modes dynamics. The sound waves dispersion relation 
can be directly obtained by linearization of Eq. (fTO]) 
around a background solution, Eq, describing the bulk 
motion of the photon fluid. By choosing the frequency 
of the driving field in order to operate at the maximum 
intracavity intensity (resonant condition) S — — ^^-\E \ 2 
and in the eikonal approximation we obtain 
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where po = \Ed\ 2 and £1 = f2 — K • Vo is the sound mode 
frequency in the locally comoving background frame. 



FIG. 1: a) 2D optical intensity and b) phase profile (mod- 
ulo 27r) obtained by numerical integration of Eq. (| 10p driven 
by vortex profile Ed with ro = 200/im and m = 6. c) The 
corresponding velocities: c a (r) (bold- line), vo(r) (solid-line), 
|t> r (r)| (dashed-line). Parameters as in the text. 

The dispersion relation (fTTj) coincides with @ except for 
the imaginary term iT which implies dissipation. This 
is not surprising since Eq. (|10|) is a dissipative system 
having the NSE as its conservative limit (T — > 0). As a 
consequence, Eq. (fTTj) exhibits the usual breaking of the 
Lorentz invariance at short wavelengths while the ampli- 
tude of low-energy modes drops to 1 /e of its initial value 
in correspondence of the propagation distance Ld — c s /T. 
However, excitations of wavelengths £ <C A <C Ld behave 
as genuine sound waves propagating with constant group 
velocity c s in direct analogy to scalar fields in curved 
spacetime. In this range indeed Eq. (fTTj) takes the form 
g^K^Ky = 0, where = (f2/c s ,K) is the covariant 
wave vector and g^ v = (g^)^ 1 - 



Experimental Proposal 

The first step of an experimental investigation would 
be to study how the presence of a "black- hole" back- 
ground modifies the propagation and the wavelength of 
sound waves. In particular, inside the the horizon, the 
background flow speed is larger than the local speed of 
sound, and so sound waves will be dragged inwards. In 
order to perform this experiment we need to apply two 
incident coherent beams to the nonlinear cavity: a first 
beam resonant with the cavity forms the photon fluid 
and creates a suitable black-hole background, and a sec- 
ond narrow amplitude-modulated beam which is modu- 
lated at the desired sound wave frequency provides a lo- 
cal excitation in a given point of the transverse plane. 
This perturbation will induce sound-like waves in the 
photon fluid which propagate away from the point of in- 
jection in the (x, y)-plane. Such waves can be detected 
by imaging the output mirror face and simultaneously 
monitoring the intensity signal at different positions of 
the transverse plane by fast photodetectors. The cor- 
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responding phase dynamics can be reconstructed by in- 
terference with a reference beam. The velocity and the 
dispersion relation of sound-waves can be obtained from 
this kind of measurements and the validity of the rela- 
tion pip can be experimentally tested. In order to simu- 
late a spacetime geometry analogue to that of a rotating 
black hole (Kerr metric), we choose the driving field pro- 
file Ed — y/pd~exp(im9 — 2iir^J^), where pd is constant. 

While numerous ex per imental methods are available to 
obtain these fields [25|, the most convenient approach 
is probably to use a phase-only spatial light modulator. 
This device consist of a two-dimensional array of indi- 
vidually addressable pixels, acting as electrically control- 
lable wave plates. A computer controls the voltage ap- 
plied to each pixel allowing the possibility to encode a 
given phase pattern into the laser beam. As an example, 
consider a 0.1-m-long cavity with T = 2 1CP 2 filled with 
85 Rb vapor at 80°C corresponding to 10 12 atoms/cm 3 . 
The self-defocusing regime is obtained by detuning the 
laser to the red side of the hyperfine transition 5Si/ 2 {F = 
2) - 5P 3/2 (F = 3) (A - 780 nm) [3, [3]. The detun- 
ing must be chosen in order to have the strongest non- 
linearity compatible with an absorption lower than T. 
A refractive index change An = n^po = 10~ 6 leads to 
a sound speed c 8 = 3 10 5 m/s. For these parameters 
£ = 390/im and Ld = 10mm, hence phononic modes are 
virtually undamped over the beam transverse dimensions 
(that can be chosen to be ss 4 mm). Fixed this value of 
c s , a sound- like wave with A = 1mm can be generated 
by a 300 MHz-amplitudc modulation of the perturbation 
beam, which can be easily achieved by an electro-optic 
modulator. Numerical integration of Eq. (fT0|) with the 
driving field £d shows that, after a short transient, a sta- 
tionary stable vortex profile forms. The intracavity field 
has a spiral equiphase surface and at the vortex centre, 
where the phase is singular, the intensity vanishes (see 
Fig. Qli,b). The corresponding sound speed and flow ve- 
locity profiles reported in Fig. [TJ: define the locations of 
the ergosurface and event horizon. The former is found 
where the flow goes supersonic, i.e. te ~ 1.12mm, while 
the latter appears where the sound speed equals the ra- 
dial component of the fluid velocity, i.e. ru ~ 0.84mm. 



A closed-form analytical expression for the intracavity 
field does not exist; however, far from the vortex core, 
where the long-wavelength approximation is valid, Eq 
asymptotes to the homogeneous solution E = In 
this case the acoustic line element is given by (8) with 

and ve 



the outer region (£ 2 /r$ < r < oo) 
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location of the horizon and ergosurface can be explicitly 
estimated. The similarity with the equatorial slice of the 
Kerr geometry becomes clearer through the transforma- 
tions of the time and the azimuthal angle coordinates in 
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After a rescaling of the time coordinate by c s the metric 
takes the form 



ds 2 oc -[l-^-(^) 2 }dt 2 +(l-—)- 1 dr 2 -2^d6dt+(rd6y 
rgr Trr ror n 



As in the Kerr spacetime the metric has a coordinate 
singularity where the radial component g rr goes to infin- 
ity which corresponds to the event horizon, i.e. rjj = 
£ 2 /r . The radius of the ergosphere is given by the 
vanishing of the temporal component of the metric gtt, 
rE = \[ r H + \Aij + 4rn 2 r#-ro/7r 2 ]. Within this surface 
gu is negative, i.e acts like a purely spatial metric coef- 
ficient, so that no state of rest can be defined. For the 
values of £ and ro used before we obtain rE ~ 1.22mm 
and rjj ~ 0.76mm in good agreement with numerical 
calculations. 



DISCUSSION AND FUTURE PERSPECTIVES 

In view of analog-gravity experiments the 2D pho- 
ton fluid presents important advantages respect to BEC 
flows. In the latters, a stationary vortex with a sink is 
not easily achieved experimentally as the matter must 
be continuously coupled out from the vortex origin, thus 
requiring some means of replenishment or the use of a 
very large condensate. On the other hand, a nonlin- 
ear optical cavity naturally support such configuration 
in the presence of a suitable driving field profile. The 
resulting steady optical vortex flow possesses two cru- 
cial properties of rotating black holes, i.e. an ergoregion 
and an event horizon. These features enables laboratory 
tests of intriguing phenomena, such as superradiance and 
Hawking radiation, which are insensitive to whether or 
not the metric satisfies the Einstein equation and robust 
against the high-energy breaking of Lorentz invariancc 
|26j |. The search of such effects will be simplified since 
light signals can be amplified and detected much easier 
than acoustic disturbances in a real fluid. In particu- 
lar, the Hawking process has a classical interpretation in 
terms of the power spectrum of exponentially redshifted 
waves emitted near the horizon which obeys the Planck 
distribution [531 . The conversion of positive- frequency 
modes to negative-frequency modes which is the basis of 
this effect has been recently observed (28|. Since in the 
photon fluid sound waves can be locally excited by inject- 
ing a second narrow beam into the cavity and detected 
by fast photodetectors, their power spectrum in the pres- 
ence of an event horizon can be measured. The observa- 
tion of a Planck distribution, although has nothing to do 
with quantum particle creation, would represent a first 
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step towards the simulation of Hawking radiation in the 
laboratory. Regarding quantum Hawking radiation, the 
possibility to detect it in such a system requires further 
investigation. 

In conclusion, sound-like waves propagation in a 2D 
photon fluid has been studied as an analog for field prop- 
agation in a curved spacetime. Although the simple light- 
matter interaction in self-defocusing media gives rise to 
an effective curved spacetime, the cavity configuration 
provides a way to control its geometry allowing the for- 
mation of acoustic ergoregions and event horizons. The 
possibility to observe phenomena analogous to superradi- 
ance and Hawking radiation and the design of a feasible 
experimental setup are currently under investigation. 
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